Absorbing state phase transitions with a non-accessible vacuum 
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We analyze from the renormalization group perspective a universality class of reaction-diffusion 
systems with absorbing states. In this class, models where the vacuum state is not accessible are 
represented as the set of reactions 2A — + A together with creation processes of the form A — ► nA 
with n > 2. This class includes the (exactly solvable in one-dimension) reversible model 2A ^iasa 
particular example, as well as many other non-reversible sets of reactions, proving that reversibility 
is not the main feature of this class as previously thought. By using field theoretical techniques we 
show that the critical point appears at zero creation-rate (in accordance with known results for the 
reversible case) and it is controlled by the well known pair-coagulation renormalization group fixed 
point, with non-trivial exactly computable critical exponents in any dimension. Finally, we report 
on Monte-Carlo simulations, confirming the field theoretical predictions in one and two dimensions 
for various reversible and non-reversible sets of reactions. 
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INTRODUCTION 

In a recent paper Elgart and Kamenev [l| have pro- 
posed a classification of absorbing state phase transitions, 
a subject that has been one of the central pillars of non- 
equilibrium statistical mechanics over the last decade 
The strategy they follow is elegant and powerful. 
The main idea is to (i) write down using standard tech- 
niques the generating functional (or, equivalently, the ef- 
fective Hamiltonian) for a given reaction-diffusion sys- 
tem; (ii) inspect the phase space in saddle-point approx- 
imation paying special attention to the "zero-energy" 
manifolds which determine the topological properties; 
(iii) detect possible structural changes in the phase por- 
trait: the birthmark of phase transitions, and (iv) classify 
them according to basic topological properties. This pro- 
cedure is a natural extension to non-equilibrium problems 
of the rearrangement of thermodynamic-potential min- 
ima occurring at equilibrium phase transitions. Hence, it 
allows for a categorization of universality classes attend- 
ing to symmetry principles, conservation laws, and few 
other relevant ingredients, which determine the phase- 
space topology and its possible structural changes. Es- 
tablishing the limits of validity of the saddle-point ap- 
proximation within this context and developing system- 
atic improvements to it remain as fundamental open 
problems. 

Using this strategy, Elgart and Kamenev report on 5 
non-trivial universality classes with absorbing states, oc- 
curring in one-dimensional systems with just one type of 
particle 0, la]- The first 4 ones are: (i) directed percolation 
(DP) characterizing generic systems with an absorbing 
phase transition and without extra symmetries, conserva- 
tion laws, quenched disorder, nor long-range interactions 
[S IB], (h) the usually called parity conserved (PC) Q 



also known as DP 2 or generalized voter class [8( which 
includes two symmetric absorbing states, (iii) the very 
elusive pair-contact-process-with diffusion (PCPD) class 
in which all reactions involve pairs of particles 0, Ho] , and 
(iv) the triplet- contact-process-with diffusion (TCPD) in 
which reactions involve triplets of particles [111]. 

In this paper we focus on the fifth class in [l|. It de- 
scribes the reversible reactions A — > 2A and 2A — > A 
occurring at rates \x and a respectively. This model was 
solved exactly in one dimension more than twenty years 
ago in a seminal paper by Burschka, Doering, and ben- 
Avraham [l2j ] by employing the empty interval method 
(l3| . Finite-size properties, scaling functions, and criti- 
cal exponents have also been exactly computed for this 
reversible model and for variations of it [121 ] in one di- 
mension. Note that except for the absence of one-particle 
spontaneous annihilation, A — > 0, this set of reactions co- 
incides with the contact process [2J a well-known model 
in the robust DP class [3j. It is, therefore, interesting 
to elucidate which is the main relevant difference in the 
renormalization group sense, giving rise to a non-trivial 
non-DP type of scaling. From considerations in [lj it 
seems that the fact that the reactions are reversible plays 
such a relevant role, but as we will illustrate, reversibility 
is a sufficient, but not a necessary, condition. 

From the field theoretical point of view, Cardy and 
Tauber had obtained in their seminal article [bj] a one- 
loop calculation of critical exponents for the closely re- 
lated set of reactions 2 A «-> 0, while in a recent paper 
Jack, Mayer, and Sollich have shown that such one-loop 
results are also valid for 2A «-> A and have to be exact ow- 
ing to the existence of detailed balance for reversible reac- 
tions [l5j . Therefore, two or more loop corrections should 
cancel out, even if this is not explicitely shown in |l5j ]. In 
any case, the main results are that the critical point is lo- 
cated at /i c = (any non-vanishing branching rate leads 
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to sustained activity) and the order-parameter critical 
exponent is (3 = 1. The long-time long-distance proper- 
ties turn out to be controlled by the well-known "pure" 
pair-coag ulation (2 A -» A) RG fixed point 0, H El E3 
and all exponents can be computed in any dimension. 

In this paper, we perform a full diagrammatic expan- 
sion of various reaction-diffusion models extending pre- 
vious analyses to all orders in perturbation theory. First, 
we recover the previously known results for the reversible 
model 2A «-> A. Afterward, using the intuition devel- 
oped from the previous full diagrammatic analysis we 
construct different sets of non-reversible reactions, and 
argue that they belong to this same universality class. Its 
key ingredient turns out to be the absence of an accessible 
vacuum state, i.e. there is no reaction mA — > but just 
pair-coagulation, combined with creation reactions of the 
form A — ► nA. The reversible reaction, n — 2 discussed in 
Q and [H is just a representative of this broader class: 
reversibility (wh ich tantamount to the detailed-balance 
condition in [151 ]) is a sufficient but not a necessary re- 
quirement. 

Let us remark that reactions as 2A «-» and its non- 
reversible extensions 2A — + 0, — > 3A, 4A, ... can also be 
argued to belong to this same class. In these cases, the 
vacuum state is accessible, but it is not stable, so they 
are not genuine absorbing state models. 

To verify the field-theoretical predictions we perform 
Monte-Carlo simulations for various non-reversible sets 
of reactions, implemented with and without hard-core 
exclusion ( "fermionic" or "bosonic" , respectively) in one 
and two dimensions. All critical exponents, are in perfect 
agreement with the RG predictions, confirming the exis- 
tence of a robust universality class, broader than thought 
before. 

Before proceeding, we should underline that while 
many of the results contained in this paper are already 
known (some from exact solutions of the reversible model 
in one dimension [lij and some from similar pertur- 
bative calculations combined with symmetry arguments 
P, [3, EH, Ell, [13]) , a systematic presentation of them, 
focusing the attention on universality aspects is, to the 
best of our knowledge, lacking in the literature. This pa- 
per aims at filling this empty space and at providing a 
comprehensive picture of this universality class, extend- 
ing it to non-reversible reactions without an accessible 
vacuum state. 



FIELD THEORY ANALYSIS OF 2A <-> A 

The techniques employed in this section are standard 
and we refer the reader to 18, 1|| and more specifically 
to 0, EE E3] for more detailed calculations and/or ped- 
agogical presentations. 

Let us apply the Doi-Peliti formalism 14, [13, [ll| (see 



and 2A — > A occurring at rates fi and a respectively. 
They can be cast into a generating functional whose as- 
sociated (bosonic) action is 

S[<t>,ir] = J dt J d d x[n(d t (t)- DV 2 (f)) - H[<f>,ir]], with 

H[<P,k] = (7T 2 -7r)(/i0-^ 2 ), (1) 

where 0(x, t) and 7r(x, t) are the density and the response 
fields respectively (some spatial and time dependences 
have been omitted for simplicity). For a general process 
kA — > jA with k and j integer numbers, the associated 
effective Hamiltonian in this formalism includes a factor 
[tt^ — tt ](/) , which is proportional to [tt 2 — tt] if and only 
if the absorbing state is not accessible, i.e. j > and 
k > 0. 

For readers with more intuition in terms of stochastic 
equations, an associated Langevin equation can be easily 
derived: 



also 2(J 2l[) to the reversible set of reactions A — > 2 A 



3t0(x, t) = DV 2 <f) + fJ,(p- ac/) 2 + y/^ - cr</> 2 r/(x, t) (2) 

where rj(x, t) is a Gaussian white noise. Let us emphasize 
the similarity between Eq. ^ and the Langevin equation 
for the DP class 0,0]. Despite of this alikeness, Eq.([2|) 
is not free from interpretation difficulties as the density 
field is not a real-valued one, but develops an imaginary 
part [2~0 ]. For this reason we avoid using it and center 
the forthcoming discussion on Eq.|T|). 

Owing to the fact that the effective Hamiltonian in 
Eq.([T]), H[(f>, tt], is proportional to (tt 2 — tt), tt — and 
7r = 1 are zero-energy manifolds. The existence of these 
two constant-7r solutions is, according to [H, at the basis 
of the non-DP behavior of this model. Indeed, the four 
zero-energy solutions: 7r = 0, tt = 1, <f> = 0, and <f> = (x/<r, 
define a rectangular geometry in the phase portrait (see 
Fig. lb and [l(), which should be compared with the 
standard, triangular, DP topology, for which only one 
constant-7r solution exists (see Fig. la and [l[) as we 
illustrate now. 

It is worth noticing that the common factor (tt 2 — tt) 
in Eq.([T]), arising from the fact that the absorbing state 
is not accessible, can be interpreted as a subtle symme- 
try between all the coefficients of (noise) terms propor- 
tional to tt 2 and their corresponding (deterministic) ones 
proportional to — tt. Indeed, it is closely related to the 
detailed-balance symmetry discussed in fl5j . If an ad- 
ditional reaction A — > occurring at rate A is switched 
on, a term A0(1 — tt) has to be added to the Hamilto- 
nian. In such a case, (tt 2 — tt) is not a common factor, 
the subtle symmetry is broken and tt = is not a zero- 
energy solution anymore. This leads to the triangular 
topology for zero-energy manifolds in the phase space 
(Fig la) and, hence, to DP-scaling. Something similar 
occurs by switching on any other reaction as mA — > 0, 
with to > 2, converting the vacuum into an accessible 
state. 
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FIG. 1: Schematic zero-energy manifolds (bold lines) for dif- 
ferent reactions. The line marked in red (non-trivial manifold, 
depending on the control parameter fi) moves downward upon 
approaching the critical point in all cases: a = j-i/a. While 
directed percolation is characterized by a triangular structure 
as in (a), models without an accessible vacuum state have 
a different form, being rectangular (b), trapezoidal (c), or 
more complicated geometries (d), for different reactions. For 
reversible models with an unstable vacuum, cj) = is not a 
zero-energy manifold (e). The structures in (c) and (d) are 
not robust under RG flow, but evolve to non-closed topologies 
as the one in figure (f). 



Let us present a different argument leading to the same 
conclusion. From standard naive power counting and rel- 
evance arguments one could be tempted to conclude that 
this problem is in the DP class, and that the critical di- 
mension is d c = 4. Indeed, as said before, the leading 
(lowest order) terms in both the deterministic part and 
the noise are identical to their analogous ones in the DP 
field theory 0,01 , for which d c = 4. The only way out of 
this naive (and wrong) conclusion is that, at the critical 
point where the linear-deterministic term coefficient van- 
ishes, the coefficient of the leading DP-like noise term tt 2 4> 
also vanishes owing to the abovementioned subtle sym- 
metry. This opens the door for higher order noise terms 
to control the (non-DP) scaling. Indeed, a proper power 
counting analysis reveals that, as the interaction Hamil- 
tonian is proportional to (tt 2 — w) , it has to be dimension- 
less, which leads to [</>] = A d (where A has dimensions of 
momentum) to ensure a dimensionless action, and con- 
sequently to [/j] = A 2 and [a] = A 2 ~ d . Therefore, the 
theory upper critical dimension is d c = 2 1, 

The existence of the common factor (tt 2 — tt) in Eq. |T]) 
implies that the /x-dependent non-trivial manifold and 
the trivial one, <j> — 0, merge at the critical point rather 
than intersecting in just a point as in DP. This is the key 



reason for the models without an accessible vacuum to 
exhibit a different type of scaling. To substantiate this 
assertion we need to prove that the previous bare-action 
symmetry, or associated topological structure, survives to 
the inclusion of fluctuations (i.e. it remains valid beyond 
mean- field [22|). 

The basic elements for a complete perturbative expan- 
sion at a diagrammatic level are: the usual DP-noise 
vertex (ir 2 (j)) [g, 3, T3] , the pair-coagulation ones (tt<P 2 
and k 2 4) 2 ) , as well as the propagator (— iu + Dn 2 + fi) 1 
l4 llrj . 1 1 71 ] . Diagrammatically, we represent response 



fields by wavy lines and the density fields by straight 
ones. For instance, the pair-coagulation noise vertex is 

depicted as a 2 ^yC^ an< ^ ana l°g° us figures, with 

different numbers of straight and/or wavy lines are em- 
ployed for the other vertices and the propagator. 

In order to perform a sound perturbative expansion 
to all orders in perturbation theory, we choose to write 
separately diagrams with and without corrections pro- 
portional to /x. The second group includes only dia- 
grams with vertices proportional to a (i.e. expansions 
of Z-functions in powers of a) . Simple inspection reveals 
that such diagrammatic corrections are those of the pair- 
coagulation process, a theory well-known to be super- 
renormalizable, i.e., all these diagrams can be computed 
and resummed to all orders. Indeed, the only possible di- 
agrammatic corrections to the pair-coagulation vertices, 
proportional to a, have the typical "bubble" structure, 
leading to a geometric series 
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l + cr 



where an is the renormalized (or "dressed" ) coagula- 
tion coefficient. Omitting external legs in Eq.Q, ctr = 
yq^y = er(l — crS) with X = x + aI , where I denotes the 
one-loop diagram evaluated at zero external frequency 
and arbitrary momentum scale A 14, 18, lj|: 
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(2tt)(* 
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A" 



Dk 2 -iui + DK 2 



(3) 

with e = 2 — d. Similar expressions are obtained for all 
the renormalized coefficients by just changing the left- 
most and/or the rightmost vertex of the series in Eq.([3]). 
If these were the only corrections, (i.e. if the RG fixed 
point was at /x c = so that diagrams including correc- 
tions proportional to /i would not give any non- vanishing 
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contribution) then the renormalized parameters would 
be: 



Mfl = = M(l - cE) 

a R = Z a o = it(1 - ctE). 



(4) 



for the two coefficients proportional to /j, and the two 
proportional to a respectively, showing that the subtle 
symmetry is not broken. The corresponding flow equa- 
tions would be 



di^R = M( 2 - cr<9;E) 
dicr R = a(e - adiT) 



(5) 



where di stands for the logarithmic derivative with re- 
spect to the momentum scale at which integral are eval- 
uated [H El 0, El- For e < 0, i.e. d > 2, the triv- 
ial (mean field) solution a = fi = is infrared stable, 
while for e > 0, the only infrared stable fixed point is 
cr* = e/(djE) with fi = 0. Plugging this into the first 
equation in ([5]) , we obtain the anomalous scaling dimen- 
sion of /ir, [i-Ir] = 2 — e = d, which coincides with the 
one-loop result obtained in [3, EH (see also [3]). 

The change in the scaling dimension of the "mass" 
term, from its naive value [fi] — 2 to the renormalized 
exact one, \\xr\ = d, induces a change in all critical 
exponents corresponding to magnitudes measured away 
from the critical point with respect to their correspond- 
ing mean field values. Moreover, as happens in pair- 
coagulation, there is no further renormalization required 



2311 and 
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for the fields nor the diffusion constant 
therefore all exponents can be exactly computed at any 
dimension. For instance, the scaling dimension of the 
field is A d and, hence, scales as fiR, implying (3 = 1 in 
any dimension. Using the same logic one obtains =2, 
v± = 1 for the correlation time and correlation length 
exponents, while right at the critical point z = 2. Using 
standard scaling relations, the density of particles as a 
function of time decays in one dimension with an expo- 
nent 9 — 1/2, while in d = 2 a similar calculation leads 
to logarithmic corrections and, in particular, to a decay 
\n(t)/t, while (3 remains equal to 1 



14] 



In order to prove that the fixed point with /i c = is not 
just a solution but also the only one one should consider 
all the possible diagrams (even if this can be done only 
in a symbolic form 24]), write down the 4 Z-functions, 
analogous to eq.((4]) for the 4 vertex functions (2 propor- 
tional to fi and 2 to cr in their bare form). Doing this, 
it is straightforward to check that 3 different and inde- 
pendent flow equations are obtained. The fourth one is 
not independent owing to the usual duality symmetry Q 
but this is not important for the argumentation here. As 
there are only 2 independent bare parameters, there is 
no way to find a fixed point for this set of 3 indepen- 
dent equations except for the trivial one ju, c = 0, which 
simultaneously satisfies in a trivial way the first 2 equa- 
tions, and leads back to the preceeding calculation, to 



the symmetry preserving Eq.(|4]), and to the same set of 
exponents. 

Note that in models in the DP class, where the naive 
power counting is different, with d c = 4, only 2 inde- 
pendent parameters in the flow equations need to be fine 
tuned to zero. The third one (corresponding to the high- 
est order noise coefficient) is irrelevant (flows to zero) al- 
ready at mean field level and, therefore, does not require 
fine tuning to vanish asymptotically. Hence, contrarily 
to the previous case, a non-trivial solution, fi c ^ exists 
leading to a DP fixed point. 

As pointed out in [l5|, the reversible reaction stud- 
ied here and 2A <-» share the same type of critical 
behavior. Indeed, the Hamiltonian in this latter case 
is (tt 2 — — o-<p 2 ) where, as before, /i and a are the 
creation and annihilation rate respectively. The zero- 
energy manifolds are: tt = ±1 and <fi = ±a/ fi/a (fig. 
le). They define a quadrangular structure, as the one 
described above, but in this case <p = is not an in- 
variant manifold: the vacuum state is accessible but it is 
not stable, so it is not properly an absorbing state phase 
transition. A perturbative analysis analogous to the one 
above can be done for the present case (indeed this is the 
model studied in [lj, ll5|) and leads to the same set of 
critical exponents; here the common factor 7r 2 — 1 plays 
the role of the subtle symmetry above. 

Finally, for reversible coagulation reactions involving 
triplets instead of pairs, A <-> 3 A, we obtain similar re- 
sults: vanishing critical point and exactly computable 
exponents, but the critical dimension is d c < 1 in this 
case [l]. 



EXTENSION TO NON-REVERSIBLE 
REACTIONS 

A careful but simple inspection of the arguments in 
the preceeding section leads to the conclusion that none 
of the reported results depends on the fact that the cre- 
ation reaction is of the form A — > 2A. As will be ar- 
gued in this section, most of them apply to more gen- 
eral non-reversible processes with creation reactions as 
A — » nA. For these, the creation part of the Hamiltonian 
is /u(7r n — tt)4>, which together with the pair-coagulation 
terms ct(7t 2 — 7r)0 2 guarantees that (tt 2 — tt) can be 
extracted as a common factor for non-reversible bare 
Hamiltonians, and hence tt is dimensionless and tt = 
and tt = 1 are zero-energy solutions as in the n = 2 case. 
For example for n = 3, H — [i(tt 3 — tt)<P — o(tt 2 — tt)4> 2 — 
(tt 2 — tt)\h(tt + l)(f> — o~ (f) 2 ] . The existence of such a com- 
mon factor in the bare Hamiltonian is, as explained be- 
fore, guaranteed if and only if the vacuum state is not 
accessible. 

For the family of non-reversible models with n > 2, the 
geometry of the zero-energy manifolds of the bare Hamil- 
tonian is not a rectangular one as occurs for the reversible 
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set of reactions with n = 2 For instance, for n = 3 
one obtains a trapezoidal geometry (zero-energy solutions: 
7r = 0, 7r = 1, = and = (7r + (see figure 

lc), but the overall topology is not changed. Indeed, as 
the critical point is approached the difference between 
the rectangle and the trapezium becomes negligible, and 
at criticality this manifold merges with the <f> — one. 
Analogously, for n — 4 one obtains a quadrangle with 3 



straight lines and a curved one (</> = (n 



l) M /a) 



(see fig. Id), which also becomes closer and closer to 
the horizontal line upon approaching the critical point. 
In all cases, the non-trivial /x-dependent manifold merges 
with the absorbing-state one <fi = at the critical point, 
and this constitutes the main trait of this class as will be 
illustrated here: in DP they intersect at criticality at a 
single point, in PC they intersect in one point in the ac- 
tive phase and merge at criticality [![, while in the class 
under scrutiny, they do not intersect in the active phase 
and merge at the critical point. 

Note that, as tt is dimensionless, all the different pro- 
cesses for different values of n are equally relevant at 
mean- field level (they just differ in powers of tt). As a 
consequence, the naive scaling dimensions for any n > 2 
are as in the preceeding section, leading to d c = 2. It is 
also important to realize that higher-order processes gen- 
erate effectively lower-order ones (in particular, A — > 2A 
is always generated) and all of them share the same de- 
gree of naive relevancy. The generation of lower-order 
processes induces changes in the zero-energy manifolds, 
and leads to combinations of the previous "pure" topolo- 
gies obtained for creation processes involving only one 
value of n. In order to render the theory renormalizablc, 
lower order terms have to be included in the bare Hamil- 
tonian, with coefficients proportional to /i (as they have 
to vanish as /i — ■> 0) that we call fi n . Indeed, from now on 
we study physical processes where various types of cre- 
ation events with different values of n are simultaneously 
present (in particular n = 2 is always generated). 

At a perturbative level, one can proceed as before, and 
separate corrections proportional and not proportional to 
\x. The first notorious difference with the reversible case 
is that upon renormalizing, the shape of some zero-energy 
manifolds is deformed if terms with n > 3 are present. 
Indeed, owing to the fact that the coefficients of 7r n , with 
n > 3, in these generalized processes renormalize as 



n(n — 1) 



(6) 



up to one loop [25| while the corresponding "mass" co- 
efficient renormalizes as in Eq.Q, different corrections 
are generated for these two coefficients equal at a bare 
level (therefore, the need to use different names, \i' n R and 
A'n.i?, = for the two of them, as a generalization of 
the single equation for lib, in Eq.(TJJ). Eq.© shows that 
the scaling dimensions of the non-linear term coefficients, 



Li' n R varies with n: the lower the value of n, the more 
relevant the corresponding non-linear term. 

Proceeding as before, it is straightforward to see by 
performing a perturbative expansion around d c = 2 that 
the only way to find a solution of the RG flow equations 
at any arbitrary order in perturbation theory is by fixing 
ll = 0. For instance, considering creation reactions with 
n = 2 and n — 3, one has 3 independent parameters: 
a, Li2 and LI3 and 5 independent flow equations. Hence, 
at criticality all creation rates have to vanish, and one 
recovers the fixed point and exponents in the previous 
section, so the universality class is preserved under the 
introduction of non-reversible reactions. Note that, in 
order to extend the calculation in the preceeding section, 
it has been enough to impose that all creation terms are 
proportional to /i. This ensures that all of them vanish 
at the critical point and generate no extra diagrammatic 
correction, but they do not need to be all equal as hap- 
pens in the reversible case. 

We should also emphasize that, as said before, the 
mass-like terms associated to each n-creation process are 
all equally relevant and they all renormalize as LI2, while 
the fx' n renormalize differently for n > 3 (see Eq.©) and 
hence, tt 2 — tt is not a common factor of the renormalized 
Hamiltonian, except at the critical point \i — where 
such a subtle symmetry is restored. The common factor 
or subtle symmetry invoked all along the calculation in 
the previous section, equivalent to the existence of re- 
versibility or detailed balance, is not a necessary condi- 
tion. As a consequence, the zero-energy manifold struc- 
ture is affected: the topology shown in figures lc and Id 
is not stable under the RG flow, tt = 1 is not a zero- 
energy manifold of the renormalized Hamiltonian, and 
the phase portrait structure becomes more complicated 
(see Fig. If). 

Despite of this, we observe that from the phase-portrait 
point of view, a key ingredient, not altered upon intro- 
ducing non-reversible reactions, is the fact that the non- 
trivial fi-dependent manifold and the trivial one <f> = do 
not intersect in the active phase and merge into a degen- 
erate manifold at the critical point. Therefore, the main 
ingredient of this universality class is not the reversibility 
nor the existence of a common factor in the renormal- 
ized Hamiltonian but the way in which the non-trivial 
manifold and the trivial one merge (2t3 | . In summary, re- 
versibility is a sufficient but not a necessary requirement. 

For completeness' sake let us comment on another fam- 
ily of reactions without an accessible vacuum, including 
higher-order creation reactions as kA — » (k + n)A with 
k > 2 which exhibit a different type of scaling behavior. 
These have to be complemented with higher order anni- 
hilation reactions as jA — > I A with j > k and j > I > 1 
in order to ensure the existence of a bounded stationary 
state. For instance, taking 2A — > 3A (with rate ff) as 
a creation reaction together with 2 A — > A (rate a), we 
need another annihilation reaction, as 3A — > 2A (with 
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FIG. 2: (Color online) Results of Monte-Carlo simulations 
for 2A — > A and A —* 3A implemented in a bosonic way in 
one dimension. The decay of the order parameter at critical- 
ity (fi c — 0) is given by t~ 0,5 (main plot). This result is well 
known as at /i = this model coincides with pair-coagulation. 
The order-parameter critical exponent is perfectly fitted by 
/3 = 1.00 (inset). Very similar results are obtained for the re- 
versible case, n — 2, as well as for higher-order non-reversible 
cases, as n = 4 and n = 5. 

rate A > /i) to have a well defined stationary state. For 
this case, even if tt = and 7r = 1 are constant energy 
solutions (at least at a bare level) the manifold cf> = 
is degenerated, and the non-trivial zero-energy solution, 
<j> = n/X — <7/(Att), intersects the line <j) = at a/[i and 
becomes singular at tt — 0, originating a very different 
topology from the one above. This topology corresponds 
to the PCPD class [H, 0|. Therefore, creation from pairs 
in systems without an accessible vacuum leads to a dif- 
ferent universality class. 

Finally, for non-reversible coagulation reactions involv- 
ing triplets (3A — ► A and A — ► nA) we obtain again that 
the universality class remains unchanged with respect to 
the corresponding reversible reaction (see last paragraph 
of the previous section) . 



MONTE-CARLO SIMULATIONS 

In order to verify the above field theoretical predic- 
tions we have performed Monte-Carlo simulations of the 
reversible reactions (reproducing some existing results 
[l2l ]) and, more relevantly, non-reversible set of reactions: 
2A — » A together with A — > nA with n = 3, 4, 5. We 
have considered two different implementations: a bosonic 
one in which the number of particles at every site in a 
lattice is unrestricted (which is the one directly related to 
the bosonic field theory presented here) , and a fermionic 
one with number occupancy restricted to be or 1. For 
both of them the same type of numerical experiments 
have been performed. Figure 2 shows our main results 
for n = 3 in the bosonic implementation. In the main 
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FIG. 3: (Color online) Results of Monte-Carlo simulations for 
2A — > A and A — > 3A implemented in a bosonic way in two 
dimensions. The decay of the order parameter at criticality 
(fi c = 0) is proportional to \n(t)/t (main plot). At its criti- 
cal point ^i = 0, this model coincides with pair-coagulation, 
for which this is a well known results. The order-parameter 
critical exponent is perfectly fitted by (5 — 1.00 (inset). 

body, we plot the time evolution of the order-parameter 
as a function of time for a one-dimensional lattice of size 
2 20 . A clean power-law decay is observed at fi c = with 
slope 9 = 0.500(1) in a log-log plot. This is not surprising 
as at fi — this model coincides with pair-coagulation. 
Direct measurements of the order-parameter as a func- 
tion of the distance to the critical point (upper inset) 
lead to j3 = 1.00(1). Also, from measures of the mean- 
squared distance associated with two point correlation 
functions [l5| one can easily measure z = 2 in all the 
cases under consideration. All the remaining exponents 
can be derived using standard scaling laws, providing a 
full check of the theoretical predictions for the bosonic 
model. For the fermionic model we obtain identical con- 
clusions. In d = 2 mean-field exponents with logarithmic 
corrections have been measured confirming that d c = 2 
(see figure 3). In d = 3, Jack et al. showed by means 
of Monte-Carlo simulations that the scaling is Gaussian 
as expected. 

We have also verified that for the sets of non-reversible 
reactions with n — 4 and with 5 one obtains the same 
set of critical exponents, supporting again the theoretical 
conclusions. 

SUMMARY 

We have shown using field theoretical arguments and 
verified by means of Monte-Carlo simulations that all 
reaction-diffusion processes including pair coagulation 
2A — > A and creation in the form A — * nA belong to 
the same universality class, regardless of whether the re- 
actions are reversible or not. The critical point is lo- 
cated at zero creation rate, and all critical exponents are 
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controlled by the well-known pair-coagulation renormal- 
ization group fixed point and can be exactly computed. 
These conclusions are in agreement with exactly known 
results for the reversible model in one dimension [T3 |. 
The main ingredient of this class of absorbing-state tran- 
sitions is that the vacuum state is not accessible and cre- 
ation occurs from individual particles. If creation occurs 
only from pairs then scaling is as in the PCPD class while, 
as soon as a reaction making the vacuum accessible, as, 
for example, 2A — > is switched-on, the system recovers 
standard DP scaling. There are also models in this uni- 
versality class as 2A <-> where the vacuum state is ac- 
cessible but in these cases it is not stable: — > 2A, 3A, 
so they are not properly absorbing-state transitions. 

We have shown that the topology of the zero-energy 
manifolds is very important to unveil universality classes, 
but there could be many subtleties leading to surprises. 
We hope that this work fosters new studies to clarify 
some of the still-standing problems on universality in 
non-equilibrium critical phenomena. 
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richsen for a very helpful exchange of correspondence. We 
acknowledge financial support from the Spanish MEyC- 
FEDER, project FIS2005-00791 and from Junta de An- 
dalucia as group FQM-165. 
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